Observation of Fluctuation-Dissipation-Theorem Violations in a Structural Glass 
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The fluctuation-dissipation theorem (FDT), connecting dielectric susceptibility and polarization 
noise was studied in glycerol below its glass transition temperature Tg. Weak FDT violations were 
observed after a quench from just above to just below Tg , for frequencies above the a peak. Violations 
persisted up to 10^ times the thermal equilibration time of the configurational degrees of freedom 
under study, but comparable to the average relaxation time of the material. These results suggest 
that excess energy flows from slower to faster relaxing modes. 
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When a system is in thermodynamic equihbrium, its 
hnear response functions (e.g. susceptibiUty) and its 
equihbrium fluctuations are related by the weU-known 
fluctuation-dissipation theorem (FDT) [1]. Special cases 
of the FDT include Nyquist's formula [2] relating elec- 
trical resistance to voltage noise, and Einstein's relation 
[3], which links the diffusion and friction coefficients in 
Brownian motion. 

There is no reason to expect the FDT to hold for sys- 
tems out of equilibrium, and it is indeed known to fail in 
far-from-equilibrium cases such as Barkhausen noise [4] . 
Glasses, however, are systems that are out of equilibrium 
but evolve very slowly. Much has been learned by study- 
ing how response functions age, that is, depend on the 
waiting time (tu,) elapsed since a temperature quench [5]. 

Recent theoretical work has shown that the FDT is 
violated by several mean-field spin glass models below 
some critical temperature [6], and also that this viola- 
tion can be used to define an effective temperature [7]. 
Computer simulations have also found FDT violations in 
finite-range spin glasses [8] , domain growth processes [9] , 
some models of structural glasses [10], elastic strings in 
random media [11], kinetically constrained lattice gases 
[12] and Ising models with dipolar interaction [13]. In all 
cases, violations are found when the characteristic obser- 
vation time, t, is of the same order or greater than the 
age of the system, i.e. when i ^ or Lutj^j ^ 1, where 
a; ~ 1/i is a measurement frequency. 

There is considerable interest in studying FDT vio- 
lations experimentally, in part because predicted viola- 
tions are model-dependent [14]. Measurements of FDT 
violations would thus furnish a new way to test mod- 
els of the glass transition and glassy dynamics. More- 
over, recent theoretical results [15] surprisingly indicate 
that non-equilibrium FDT violations can, at least in spin 
glasses, give information about an equilibrium order pa- 
rameter, which would be otherwise inaccessible. 

Aging in susceptibility [16] and noise [17] has been 
measured in structural glasses. However, measurements 



of the relation between noise and response thus far re- 
ported in structural [18] as well as spin [19] glasses are 
compatible with the FDT within experimental errors. 
This gap between theory and experiment may be due 
to the fact that the regime t^uj ^ 1 is very difficult to 
access and/or that actual violations are rather weak. 

In this letter we present results that show a weak vi- 
olation of the FDT in a supercooled liquid, glycerol, fol- 
lowing a quench below its glass transition temperature. 
The effective temperature (Tcff) measured via dielectric 
noise slowly relaxes towards the thermal bath tempera- 
ture T. Most surprisingly, violations are observed even 
with t^uj 3> 1. The relaxation is so slow compared to 
the thermal equilibration time of the configurational de- 
grees of freedom (CDF) probed by the experiment, that 
these results can only be understood in terms of strong 
coupling between CDF. 

The experiment is similar to the "oscillator as ther- 
mometer" gedanken experiment proposed by Cugliandolo 
et al. [7]. The idea is to connect a capacitor, whose di- 
electric is the glass under study, in parallel with a loss- 
less inductor (see fig. 1). In thermodynamic equilibrium, 
the average energy of the oscillator, C{V'^), will be ksT 
by equipartition, where {V"^) is the integrated voltage 
noise, ks Boltzmann's constant and T the temperature. 
The noise power spectral density, Sy, will be given by 
Nyquist's formula [2], 

(1) 

Here Z is the circuit impedance, L the inductance, and 
C = C + iC" = Cqe the capacitance (with e the complex 
dielectric susceptibility) . 

The strongly peaked noise spectrum produced by the 
resonant circuit (see fig. 2) provides three independent 
quantities: peak position, width and height, which relate 
directly to C", C" and T, respectively. If the FDT ceases 
to be valid, then this would be apparent as a failure of 
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FIG. 1. Experimental setup. The sample forms the dielec- 
tric of a capacitor in a resonant circuit, which is driven by the 
polarization noise of the sample. 



the noise spectrum to satisfy Nyquist's equation with T 
as the thermodynamic bath temperature. Thus, the need 
for separate measurement of noise and response is elim- 
inated with this method. We note that measurement of 
broadband noise alone [17] only provides two quantities 
(spectral density and exponent) and therefore is insuffi- 
cient to test FDT violations. 

For a non-equilibrium (aging) glass, the FDT is be- 
lieved to fail [6], and an effective temperature, Teff(a;, tw), 
can be defined [7]. It should replace T in eq. 1, control the 
direction of heat flow, and determine the average energy 
of an oscillator coupled to the sample (as in the equipar- 
tition theorem). Tes{ijJ-,tw) should approach the bath 
temperature for — > oo at any fixed frequency, but its 
behavior is otherwise system-dependent. Whether this 
scenario is completely valid or not, we can use T^a to 
quantiiy any observed FDT violations. 

The sample capacitor was fabricated by rolling two 
copper foils separated by a sheet of high porosity pa- 
per on a copper tube, to which a heater and temperature 
sensor were attached. Liquid glycerol was then absorbed 
into the paper. The capacitor was connected in parallel 
with an inductor and both placed inside a liquid nitro- 
gen cryostat at low pressure (« 100 mTorr). To minimize 
the noise due to the lossy part of the inductance, the in- 
ductor was placed at the bottom of the cryostat in good 
thermal contact with the liquid nitrogen. The voltage 
noise was amplified with a Stanford SR560 preamplifier, 
filtered with a Stanford SR640 low-pass filter, digitally 
sampled with a personal computer and Fourier-analyzed 
to obtain the voltage power spectrum as a function of 
aging time. 

Since capacitor and inductor are at diflferent tempera- 
tures, cq. 1 has to be replaced by a generalized formula 
[20], which for our circuit gives 

„ , ^ 2kB -cJ^T,s{u;o,t^)\L\^C" ^ ujUL" 
^^y"^) l-2w2(L'C"-L"C")+w1i|2|C|2' ^' 
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FIG. 2. Typical voltage power spectrum for 179.8 K 

(squares), and fit to eq. 2 (solid line). Measurement are made 
after bath temperature T has stabilized. Inset: Typical ther- 
mal history. 

where To = 77K is the inductor temperature, which now 
has a lossy part {L = L' + iL"), which is measured sepa- 
rately with a known capacitor. We disregard any possible 
CO dependence of T^e because Sv is strongly peaked, and 
its value at the resonance frequency, luq, will dominate. 
All other quantities being known, a least-squares fit of 
the power spectra to eq. 2 yields C, C" and TcS- 

The sample was always heated to a bath tempera- 
ture of 210 K (above the glycerol glass temperature) and 
then cooled to T near 180 K. A typical thermal history 
is shown in fig. 2 (there is an unavoidable temperature 
undershoot). All reported results are from data taken 
after the bath temperature stabilized at its final value. 
t-u, =Q was taken as the time when the sample falls out of 
equilibrium. This depends on cooling rate and thermal 
history. We have determined Tg to be w 196 K. Adjust- 
ment of Tg by a few degrees will shift the time axis by 
at most a few hundred seconds, without aff'ecting our re- 
sults. 

Fig. 2 shows a typical power spectrum, together with 
the fit using eq. 2. The spectrum is an average obtained 
from voltage noise sampled in a time window around a 
particular t^)- The time window increased with i^, so as 
to reduce errors for large tyj. For shorter t-u, (< 10^ s), it 
was necessary to average data from several (^ 15) runs. 
The resonance frequency is /o ~ 7.7 Hz. The errors in 
the fitting parameters were estimated by fitting simulated 
power spectra with Gaussian noise. 

Figure 3 plots C and C" vs. waiting time at T = 
179.8 K, which clearly show the aging of the dielectric 
susceptibility, as previously seen by Leheny and Nagel 
[16]. The curves can be fitted by a stretched exponential, 
C'{t^) =C'^ + AC exp{-(t^/r)^}, with r and P taken 
from ref. [16]. 

The effective temperature is plotted in fig. 4 against 
for T = 179.8 K. Two features of this graph are notewor- 
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FIG. 3. Real (squares) and imaginary (circles) capacitance 
vs. waiting time, showing aging in the dielectric susceptibility 
at T = 179.8 K. Solid lines: fits to stretched exponentials. 




FIG. 4. Effective temperature T^s vs. waiting time for 
T — 179.8 K. Points: experimental data, solid line: fit to 
exponential decay. 



thy. First, a fit of the equilibrium data (long tw) with 
eq. 2 yields the bath temperature with high accuracy. 
Second, a fit of eq. 2 to the non-equilibrium data at early 
times fails to yield the bath temperature. This failure to 
fit Nyquist's formula can only be explained by an FDT 
violation. Although the violations are small, the statis- 
tical certainty of a violation is very high: we performed 
the fitting procedure on simulated spectra with the same 
susceptibility as the experimental spectra but obeying 
the FDT. With T = 179.8 K, the T^s obtained by fitting 
500 000 simulated spectra was in no case greater than 
183.67 K, while 5 experimental points at the earliest 
are above this value. Tes relaxes slowly towards T on a 
time scale similar to that of the aging in susceptibility. 
The decay can be fitted equally well by an exponential 
or a stretched exponential. Most striking, this FDT vio- 
lation persists much longer than expected, to Lot^ > 10^. 
This effect was clearly observed only in a narrow range of 
temperatures around 180 K (fig. 5), while susceptibility 
aging was seen over the range 177K-184K, in agreement 
with [16]. 

It is interesting to compare and contrast the new re- 
sults for Tofr with other thermodynamic measurements. 
For example, dynamic heat capacity measurements in 
glycerol [21] show that changes in temperature are ac- 
companied by heat flow on long time scales similar to 
those found in other dynamical measurements. There 
are two ways this could be relevant to the present exper- 
iments. 

First, heat flow will lead to a temperature gradient 
across the sample, and thereby an apparent rise in T^s. 
However, if this was important, the susceptibility at early 
times should correspond to that of higher temperatures 
in equilibrium, but we find that C"(179.8K) is at all ty, 
lower than the final value of C"(181.7K). Furthermore, a 
generous estimate of the maximum temperature rise gives 
50 mK at 180 K, negligible compared with the measured 



rise in Teg. 

Second, the dynamic heat capacity measurements can 
be understood [22] in terms of an excess enthalpy, AH, 
trapped in the sample for times comparable to the re- 
laxation time. A fictivc temperature, T/, [23] can be 
defined as the temperature at which the glass and liquid 
lines intersect in a diagram of enthalpy vs. temperature. 
Then Tf = T + AH/ AC, where AC = Cg - Ci is the 
heat capacity difference between the liquid and the glass 
[22]. As the glass relaxes, its enthalpy tends to that of 
the equilibrium liquid, and so Ty T . It is tempting 
to attribute the excess Teff observed simply to the excess 
enthalpy per CDF. We can use T/ — T as a rough esti- 
mate of the excess enthalpy per CDF stored in the glass 
as a function of tw We computed T/ as described in [22] 
using the particular thermal history of our sample. As 
fig. 5 shows, the deviation of Tcff and Ty from T is similar 
in magnitude at 180 K, but the temperature dependences 
of the two are distinctly different. 

That Toff and Ty behave differently is not unexpected, 
since the oscillator only couples to those CDF with re- 
laxation time r~l/a;~0.1s. In a parallel relax- 
ation picture, excess enthalpy trapped in these relatively 
fast CDF would relax in ~ 0.1s and independently of 
slower CDF, too fast for this experiment to detect. It is 
possible that step-like relaxation of slow CDF near the 
a peak could contribute enough to be observed at the 
much higher measurement frequencies, as in Barkhausen 
noise [4]. Barkhausen noise decreases with frequency 
as 1/w^, whereas the equilibrium noise decreases more 
slowly. This would require a highly improbable T^a of 
10^ K at the a peak frequency to increase Toff by 4 K at 
the measurement frequency. Moreover, Monte Carlo sim- 
ulations [24] of a parallel kinetics model with the same 
relaxation times as our sample show no significant in- 
crease in the noise power at the measurement frequency 
during a quench. Thus our results are inconsistent with 
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FIG. 5. ATeff = Teff (t,„ = 10000 s) - %s{tn, = oo) vs. bath 
temperature (points). Solid line is Tf — T. 

parallel kinetics. 

Previous results, such as the aging of dielectric re- 
sponse [16] and dielectric hole burning [25] could be 
qualitatively understood in terms of a slow relaxation 
of relaxation rates via coupling between barriers and lo- 
cal density. But the present results further require that 
the CDF under observation be continually 'heated' above 
the phonon temperature. Hence some mechanism must 
transfer excess energy stored in slower CDF to faster 
CDF active in the experimental window. This could hap- 
pen through reshuffling of relaxation times of individual 
CDF, or by directed energy exchange between different 
CDF. The peak in figure 5 then arises naturally, because 
at high T the relaxation is too fast to be seen, while at 
low T the energy release from relaxing CDF is too slow. 
Some form of series kinetics [26] (as seen in recent sim- 
ulations [27]), in which relaxation occurs in a series of 
constrained steps, may explain these results. 

In conclusion, we have reported experimental observa- 
tion of weak FDT violations in glycerol, which we have 
analyzed in terms of an effective temperature. Like other 
quantities, the effective temperature ages, but provides 
different information than fictive temperature, suscepti- 
bility, or dynamic heat capacity. Specifically, the present 
results indicate that successful models of the glass tran- 
sition must explicitly account for nontrivial energy ex- 
change between slower and faster configurational degrees 
of freedom. Extending these results to a broader range 
of temperature, frequency, and other materials would be 
useful. 

We thank L. F. Cugliandolo and M. B. Weissman for 
helpful discussions. This work was supported by NSF 

DMR 94-58008 and 98-77069, and ACS-PRF 33329- AC5. 
TSG was partly supported by Fundacion Antorchas (Ar- 
gentina). 



* e-mail: t-grigera@usa.net 
t e-mail: israeloff@neu.edu 
[1] R. Kubo, Rep. Progr. Phys. 29, 255 (1966); R. Kubo, 

M. Toda and N. Hashitsume, Statistical Physics II (2nd 

ed.), Springer- Verlag, Berlin (1991) 
[2] H. Nyquist, Phys. Rev. 32, 110 (1928) 
[3] A. Einstein, Ann. d. Phys. 17, 549 (1905) 
[4] H. Barkhauscn, Z. Phys. 20, 401 (1919) 
[5] L. C. E. Struik, Physical ageing in amorphous polymers 

and other materials, Elsevier, Amsterdam (1978) 
[6] L. F. Cugliandolo and ,J. Kurchan, Phys. Rev. Lett. 71, 

173 (1993); J. Phys. A 27, 5749 (1994); Phil. Mag. B 71, 

501 (1995) 

[7] L. F. Cugliandolo, J. Kurchan and L. Peliti, Phys. Rev. 
E 55, 3898 (1997). 

[8] S. Franz and H. Ricgcr, J. Stat. Phys. 79, 749 (1995); 
E. Marinari, G. Parisi, F. Ricci-Tersenghi and J. J. Ruiz- 
Lorenzo, J. Phys. A 31, 2611 (1998) 

[9] A. Barrat, Phys. Rev. E 57, 3629 (1998) 
[10] G. Parisi, Phys. Rev. Lett. 79, 3660 (1997); Phil. Mag. 
B 77, 257 (1998); J.-L. Barrat and W. Kob, Europhys. 
Lett. 46, 637 (1999); J. Phys.: Condens. Mater 11, A247 
(1999) 

[11] H. Yoshino, Phys. Rev. Lett. 81, 1493 (1998) 
[12] M. Sellitto, cond-mat/9804168 

[13] D. A. Stariolo and S. A. Cannas, Phys. Rev. B 60, 3013 

(1999) 

[14] A review is J.-P. Bouchaud, L. F. Cugliandolo, J. Kur- 
chan and M. Mczard, in: A. P. Young (ed.), Spin Glasses 
and Random Fields, World Scientific, Singapore (1998) 

[15] S. Franz, M. Mczard, G. Parisi and L. Peliti, Phys. Rev. 
Lett. 81, 1758 (1998) 

[16] R. L. Leheny and S. R. Nagol, Phys. Rev. B 57, 5154 
(1998) 

[17] N. E. Israeloff and T. S. Grigera, Europhys. Lett. 43, 308 
(1998) 

[18] N. E. Israeloff, Phys. Rev. B 53, R11913 (1996) 

[19] M. Ocio, H. Bouchiat and P. Monod, J. Phys. Lett. (Fr.) 

46, L647 (1985); W. Reim, R. H. Koch, A. P. Malozc- 

moff, M. B. Ketchen and H. Maletta, Phys. Rev. Lett. 

57, 905 (1986); H. Bouchiat and M. Ocio, Comments 

Cond. Mat. Phys. 14, 163 (1988) 
[20] W. R. Bonnet, Electrical Noise, McGraw-Hill (1960) 
[21] N. O. Birge and S. R. NagcL Phys. Rev. Lett. 54, 2674 

(1985); N. O. Birge, Phys. Rev. B 34, 1631 (1986) 
[22] S. L. Simon and G. B. McKenna, J. Chem. Phys. 107, 

8678 (1997) 

[23] A. Q. Tool, J. Am. Ccrani. Soc. 29, 240 (1946); O. S. 

Narayanaswamy, ,]. Am. Coram. Soo. 54, 491 (1971); C. 

T. Moyniham at al. J. Am. Ceram. Soc. 59, 12, 16 (1976) 
[24] T. S. Grigera and N. E. Israeloff, to be published. 
[25] B. Schiener, R. Bohmer, A. Loidl and R. V. Chamberlin, 

Science 274, 752 (1996); B. Schiener, R. V. Chamberlin, 

G. Diezemann and R. Bohmer, J. Chem. Phys. 107, 7746 
(1997) 

[26] G. H. Fredrickson and H. A. Andersen, Phys. Rev. Lett. 
53, 1244 (1984); R. J. Palmer, D. L. Stein, E. Abrahams 
and P. W. Anderson, Phys. Rev. Lett. 53, 958 (1984) 

[27] C. Donati, J. F. Douglas, W. Kob, S. J. Plimpton, P. 

H. Poole and S. C. Glotzer, Phys. Rev. Lett. 80, 2338 
(1998) 



4 



